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This paper is concerned with the effect of a weak spanwise- variable mean-flow distortion 
on the growth of oblique instability waves in a Blasius boundary layer. The streamwise 
component of the distortion velocity initially grows linearly with increasing streamwise 
distance, reaches a maximum, and eventually decays through the action of viscosity. This 
decay occurs slowly and allows the distortion to destabilize the Blasius flow over a relatively 
large streamwise region. It is shown that even relatively weak distortions can cause certain 
oblique Rayleigh instability waves to grow much faster than the usual two-dimensional 
Tollmien-Schlichting waves that would be the dominant instability modes in the absence 
of the distortion. The oblique instability waves can then become large enough to interact 
nonlinearly within a common critical layer. It is shown that the resulting nonlinearity is 
weak and that the common amplitude of the interacting oblique waves is governed by the 
amplitude evolution equation derived in Goldstein & Choi (1989). The implications of these 
results for Klebanoff-type transition are discussed. 


1. Introduction 


Transition to turbulence in boundary layers usually begins with initially linear and non- 
interacting instability waves that grow to nonlinear amplitudes as they propagate down- 
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stream. The first nonlinear stage of evolution - which might more appropriately be re- 
ferred to as a modal interaction stage - is usually characterized by the rapid growth of 
three-dimensional disturbances due to resonant interactions between instability waves and 
between instability waves and streamwise vortices. 

This phenomena is usually studied experimentally by exciting the flow with relatively 
two-dimensional single-frequency excitation devices. The initial motion, say just down- 
stream of the excitation device, should then be periodic in time and reasonably well de- 
scribed by linear instability theory - provided, of course, the excitation levels are sufficiently 
small. The typical mean flow is relatively two-dimensional and fairly close to a Blasius pro- 
file at the low Mach numbers where most of the experiments have been carried out. The 
instability wave growth rates should then be small compared to the inverse of the mean 
boundary-layer thickness in these experiments. 

When flow visualization devices are used, the initial modal interaction stage is evi- 
denced by the appearance of A-shaped structures which can either be aligned or staggered 
in alternating rows. The aligned arrangement, which occurred in the original Klebanoff 
& Tidstrom (1959) and Klebanoff, Tidstrom &; Sargent (1962) experiments, is usually re- 
ferred to as ‘peak-valley 1 splitting. It is believed to be a complex phenomena (Kachanov 
& Levchenko 1984, §5.2) that can be explained in terms of at least three different (rela- 
tively weak) resonant-type interaction mechanisms each of which probably plays a role in 
one or more of the many experiments that have been carried out to study this phenomena 
(Kachanov, Kozlov & Levchenko 1985; Kachanov 1987; Hama & Nutant 1963; Kovasznay, 
Komoda & Vasudeva 1962; Nishioka, Asai & Iida 1979; and others). 

A resonant-type interaction involving weak streamwise vortices seems to have played 
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an important role in the original Klebanoff & Tidstrom (1959) and Klebanoff et al. (1962) 
experiments. The present paper is an attempt to provide a systematic first-principles analy- 
sis of this phenomena by using high-Reynolds-number asymptotic methods. There has been 
a tendency in the literature to separate such analyzes into wave-wave and wave-vortex in- 
teractions. This paper, which brings together a number of recent ideas in order to shed 
some light on the KlebanofF-type transition process, turns out to be a combination both of 
these approaches. 

We first consider the initial linear region, just downstream of the excitation device, 
where the instability waves are still small enough so that no significant modal interac- 
tions take place. The instability waves will grow on the relatively slow viscous time scale 
when the mean flow is two-dimensional, with the two-dimensional mode exhibiting the 
most rapid growth. However, even relatively weak spanwise-periodic mean-flow distortions 
(i.e. streamwise vortices) can cause certain oblique modes to grow on the inviscid time scale 
through a kind of resonant-interaction mechanism first considered for Gortler vortices by 
Nayfeh (1981) and later by Bennett & Hall (1988), Nayfeh & Al-Maaitah (1988) and Hall 
& Seddougui (1989). This resonant interaction allows the oblique modes to grow faster 
than the plane wave once the Reynolds number becomes sufficiently large. The streamwise 
vortices, which are generated whenever quasi-periodic cross-flow velocities exist in the flow, 
have streamwise velocity components that initially grow in proportion to the downstream 
distance. These velocity components can then become quite large before viscous and/or 
nonlinear effects cause them to saturate 1 . 

'This is related to, but somewhat different than, the algebraic growth mechanisms studied by Ellingsen 
& Palm (1975), Hultgren fc Gustavsson (1981) and Landahl (1990) 
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It is therefore appropriate to suppose that the spanwise-periodic motions are initiated 
by a steady cross flow with spanwise wavenumber, say 2/3, and a pair of equal-amplitude 
oblique instability modes with the same streamwise wavenumber and scaled frequency but 
opposite spanwise wavenumbers, say ±/3. These two modes combine to form a standing 
wave in the spanwise direction that propagates only in the direction of the free stream. 
This situation is typical of wave excitation experiments which often involve relatively long 
excitation devices oriented perpendicular to the free-stream direction. 

When the Reynolds number i6 sufficiently large, the streamwise vortices can persist over 
streamwise distances that are long enough to enable the oblique modes to reach nonlinear 
amplitudes. Since the oblique-mode growth rates turn out to be small compared to the 
spanwise wavenumbers in the present analysis, the initial nonlinear interactions are confined 
to a localized region centered around the so called ‘critical level’ where the streamwise 
component of the mean-flow velocity is equal to the common phase speed of the oblique 
modes. The flow outside the critical layer is still governed by linear dynamics and is given 
by the superposition of a Blasius flow, a spanwise-periodic mean-flow distortion and a pair 
of oblique instability modes. The common amplitude of the oblique modes is completely 
determined by the nonlinear dynamics within the critical layer. 

The critical-layer nonlinearity turns out to be weak in the sense that it enters through 
an inhomogeneous term in a higher-order problem rather than through a coefficient in the 
lowest-order or dominant-balance equation. This ultimately means that the oblique-mode 
amplitude is completely determined by a single amplitude-evolution equation. It turns out 
that this equation is the same as the one that was obtained by Goldstein & Choi (1989), who 
considered the related problem of the interaction of two oblique modes in a two-dimensional 
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shear layer. 

Since this type of interaction has already been shown to play a role in both the subhar- 
monic transition process (Goldstein & Lee 1992; Mankbadi, Wu & Lee 1993; and Wundrow, 
Hultgren & Goldstein 1994) as well as in one of the major competing scenarios for the har- 
monic transition process (Goldstein & Lee 1992), this adds to the growing evidence that it 
is a relatively universal mechanism that can occur in many of the transition processes iden- 
tified in the literature. It is our belief that this identification of universal mechanisms is one 
of the important ways in which the fundamental theory can contribute to our understanding 
of the very complex and multifaceted transition process. 

The paper is organized as follows. The spanwise-periodic mean-flow distortion is an- 
alyzed in §3. The linear stability of this flow is considered in §§4 and 5 where it is shown 
that even relatively small distortions can cause the oblique instability modes to grow more 
rapidly than the fastest growing plane wave when the Reynolds number is sufficiently large. 
The nonlinear effects are considered in §6 where it is shown that the nonlinearity is weak 
and the common amplitude of the oblique modes is governed by the amplitude evolution- 
equation derived in Goldstein & Choi (1989), but with different numeric coefficients. The 
implications of the results are discussed in §7. 

2. Formulation 

To fix ideas, we consider the incompressible flow over an infinitely thin flat plate and 
suppose that a small-amplitude spanwise-periodic motion is suddenly imposed on the flow 
(say by a vibrating ribbon or other excitation device) at a distance L downstream from the 
leading edge. The Cartesian coordinates system ( x,y,z ) is attached to the plate at L with 
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x aligned with the free stream, y normal to the plate, and 2 in the spanwise direction. All 
lengths are non-dimensionalized by ft where 

6 = L/R i (2.1) 

is characteristic of the mean boundary-layer thickness at L and 

R = LUoo/v (2.2) 

is the global Reynolds number based on the free-stream velocity Uoo and the kinematic 
viscosity v. The time t, velocity u = (u,v,w), and pressure variation p from the free- 
stream value Poo are non-dimensionalized by 6 /Uoo, Uoo and pU^, respectively, where p is 
the density. With this non-dimensionalization, the Navier-Stokes equations become 


(d/dt + u-V) u + Vp = , 

(2.3) 

0 

II 

3 

> 

(2.4) 


where V = id/dx + jd/dy + kd/dz is the gradient operator. 

When the excitation device is placed perpendicular to the free-stream direction, as 
we shall now assume, it will generate a disturbance that consists of a weak mean-flow 
distortion plus a pair of equal-amplitude oblique instability modes that form a standing 
wave in the spanwise direction. Since we ultimately require that the spanwise length scale 
of the disturbance, say 1//3, be much larger than the local boundary-layer thickness, we 
introduce the scale factor <7 through 

<tP = 6(3, 0 < a sj 1, /3 = 0(1), (2.5) 

and eventually require that a -C 1. The distinguished scaling for the imposed cross-flow 
velocity at x = 0, say oR~% Wo(y,az) where a = 1/lna, then corresponds to Wo = 0(1). 
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3. Mean-flow distortion 

The imposed cross flow dR~^Wo(y, az) generates a mean-flow distortion that ulti- 
mately decays out through the action of viscosity, but causes the mean flow to become 
inflectional before this occurs. This allows certain initially linear instabilities to grow to 
nonlinear amplitudes by essentially inviscid mechanisms. It turns out that the most inflec- 
tional profiles lie in the region where 

x 2 = x/cr 3 R$ = 0(1) (3.1) 

for 0 < a ^ 1. This scale will be very long compared to the triple-deck length scale when 

<r>l/i?i (3.2) 

which is now assumed to be the case. When a = 0(1), the mean-flow velocity U = ( U , V, W) 
and pressure P expand like 


U = (U,R-lV,R~lW) + (3.3) 

P = R~ 1 P + ---, (3.4) 

and are determined by the parabolized Navier-Stokes equations 

(Ud/dx 2 + U-V T )U + V T P = S7\V, (3.5) 

dU/dx 2 + V t -17 = 0, (3.6) 


where U = (U , V, W) and V r = jd/dy+kd/dz is the gradient in the transverse plane; but, 
when a <C 1, the mean flow will be interactive relative to the cross-flow direction and will be 
determined by the three-dimensional boundary-layer equations in the main part of the flow 
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field. The initial distortion becomes fully interactive in a sufficiently small neighborhood 
of x = 0, but this region is of little interest here since (with the present scaling) it has 
negligible effect on the instability waves. Its structure is discussed in appendix B for the 
long- wavelength limit a ■< 1 in order to show how the x 2 -scale flow (to be discussed below) 
evolves from the initial disturbance. Since (3.5) and (3.6) must be solved numerically, 
the relevant physical mechanisms can best be understood by concentrating on the long- 
wavelength limit a < 1 for which analytic solutions can be obtained. The discussion of the 
order-one-wavelength case is deferred to the end of §7. The structure of the long- wavelength 
solution is similar to the short-wavelength triple-deck solutions worked out by Rozhko & 
Ruban (1987) and by Choudhari, Hall & Streett (1992), but some new results are obtained. 

3.1. The main boundary layer 

In the main region where x 2 and y are both order one, the spanwise- variable mean flow is an 
inviscid perturbation about the local Blasius profile (U B , R~$ V B ). The mean-flow velocity 


and pressure in this region expand like 

U = U B + o 4 U d + 0(<r 5 ), (3.7) 

V = R-*[V b + ctV d + 0(ct 2 )}, (3.8) 

W = <jR~^\W d + 0(<rlncr)], (3.9) 

P = <r 3 R- l [P D + 0(o)] t ( 3 . 10 ) 

where 

U B (x 2 ,y) = f'(y) - (T 3 \x 2 yf"(y), V B (y) = \ [yf'(y) - f(y )\ , (3.11) 
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/ denotes the Blasius function which satisfies 


f" + W = 0, m = /'(O) = o, /'(oo) = 1 , (3.12) 

t/ D , Vd, and W D are functions of x 2 , y and 


z = oz, (3.13) 

P D is a function of x 2 and z only, and a prime denotes differentiation with respect to the 
argument. Substituting (3.7)-(3.10) into (2.3) and (2.4) and using the fact that 


W D -* W 0 (y, z) as x 2 -*■ 0 


leads to 


Vo = f" (A + crx 2 H s ) , 
V D = - f'(A X3 +dH;), 
W D = Wo, 


(3.14) 


(3.15) 

(3.16) 

(3.17) 


where 


rr w 0 (s,z) wq(o,z) 


H{y,z) = / 
Jo 


, , WU), 

ds + 7 In y, 


(3.18) 


L f'(s) a o5 

Aq « 0.33206 is the Blasius wall-shear stress and an independent variable used as a subscript 
denotes differentiation with respect to that variable. The particle displacement A and 
pressure distribution P D are, at this point, arbitrary functions of x 2 and z, which have the 
implicit a dependence 


A = Aq{x 2 ,z) + aAi(x 2 ,z), P D = P 0 {x 2 ,z) + aPi(x 2 ,z), 


(3.19) 
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where the 0(1) terms have been inserted to facilitate matching with the solution in the 
viscous wall layer to be discussed below. We do, however, require that 

A, Ax 2 -* 0 as x 2 — ► 0 (3.20) 

so the velocity distortions U D and V D arise from only the imposed cross flow Wo. A a 
dependence similar to that in P D will be implicitly assumed in the corresponding dependent 
variables in the wall layer as well as in the inviscid outer layer which is considered next. 

3.2. The outer layer 

Since the solution (3.16) does not vanish as y —* oo, it is necessary to introduce an outer 
region where 

y = oy = 0(l) (3.21) 

in order to satisfy the appropriate free-stream boundary conditions. The solution in this 
region relates P D and A. 

Substituting (3.1), (3.13) and (3.21) into (2.3) and (2.4) and matching with the main- 
boundary-layer flow shows that the solution in the outer layer expands like 


U = 1 + <J- Z R~ X U + ■■■, 

(3.22) 

V = friv fl (oo) + <TiriF + -.., 

(3.23) 

W = oR-*iW + •••, 

(3.24) 

P = o~*R- l P+ •••, 

(3.25) 


where it has been assumed that Wo = o(y *) as y — » oo. The functions U , V, W and 
P of X 2 , y and z are determined by the linearized Euler equations. It follows from these 
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equations that 


Pyy + Pzz = 0 (3.26) 

while the free-stream boundary conditions and matching with the main-boundary-layer 
solution require that 

P — * 0 as y -+ oo, (3.27) 

and 

P = P D , Py = Ax 2 x 2 at y = 0. (3.28) 

The solution to (3.26)-(3.28) is most easily found by expressing P as the real part of an 
analytic function of the complex variable z -f i y. The Cauchy integral formula can then be 
used to relate P D and A. For the spanwise-periodic mean flows that are of interest here, 
this leads to 

8 r n /P 

Pdz = - T A X2 x 2 {x2,(,\°)zoiP(z - C)dC, (3.29) 

where -f denotes the Cauchy principal value. 

3.3. The viscous wall layer 

It will now be assumed that 

B(z) = W o {0,z)t0. (3.30) 

(The case where Wo(0, z) = 0 is much more complex and is best studied by considering 
the order-one-wavelength problem discussed at the end of §7.) It then follows from (3.15), 
(3.18) and (3.30) that 

Ud —* + 0 X 2 B' In y as y — ► 0 (3.31) 
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and consequently that U will not satisfy the proper wall boundary condition. It is therefore 


necessary to introduce a viscous wall layer where 


II 

VS 

ii 

o 

(3.32) 

in order to bring U to zero at the wall. The mean velocity in this region expands like 

U = ct\ 0 Y + <t\U - ±\ 2 0 Y 4 - |A 0 a: 2 r) + • • • , 

(3.33) 

V = <7 2 R-i{V +\\ 0 Y 2 ) + ---, 

(3.34) 

W = R~lW + ---, 

(3.35) 

where U, V and W are functions of x 2 , Y and z. 


Substituting (3.33)-(3.35) and (3.10) into (2.3) and (2.4) yields 


A 0 YU X2 + X 0 V = U Y y , 

(3.36) 

A 0 YW Xa + P Di = Wyy, 

(3.37) 

v X2 + v Y + W t = 0 . 

(3.38) 

It follows from (3.30) and (3.31) that the solution to these equations 

must satisfy the 

boundary conditions 


U = V = W = 0 at y = 0, 

(3.39) 

and 


U Ao.4 4- crxiB' ln(<ry ), W — » bB as Y — » 00 . 

(3.40) 

The solution must also satisfy the upstream condition 


U — > 0, W — > dB as X 2 — > 0, 

(3.41) 
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since the thickness of the viscous wall layer goes to zero as 12 0. 

Eliminating V between (3.36) and (3.38) leads to 

Ao YUyx 2 ~ AoM's = Uyyy, (3.42) 

which determines U once W is known. This equation must be solved subject to 

Uyy =0 at Y = 0 (3.43) 

in addition to (3.39)-(3.41). Since the boundary-value problem (3.37) and (3.39)-(3.43) 
has a similarity solution when P D is set equal to zero and A is selected appropriately, it is 
convenient to seek a solution of the form 

U = ax 2 B'(z)F(r,) + t/ (,) (x 2 ,y, 2; a ), W = aB(z)G{r]) + W' (,) (z 2 , 5; a) (3.44) 

where 

V = (Ao /x 2 )iY, (3.45) 

and F and G are determined by 

F"' +\ri i F" -\i)F' = -G, G" + Aq 2 G' = 0 (3.46) 

with 

F = F" = G = 0 at x] = 0, (3.47) 

and 

F —* In T], G — > 1 as r) —* oo. (3.48) 

It follows that 

G = r(J)- , 7 (i, r ), (3-«) 
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where r = t? 3 /9, V is the usual gamma function and 7 is the incomplete gamma function 
defined on page 260 of Abramowitz &: Stegun (1964). It can also be shown that 


F = |r(i)r(lr l £ - fr(l)-' [r».-' + (r + 2)7(1, r)] , (3.50) 

where M denotes the confluent hypergeometric function in the notation of Abramowitz & 
Stegun {1964, p. 504). 

The inhomogeneous solutions f/ (/) and W (I) can be found by taking Laplace transforms 
with respect to z 2 . The details are given in appendix A where it is shown that 

xPd ss = -r(i)Q? } - nAx(Z,z;o)(x 2 -0"*d£, (3-51) 

JO 

and 

4xU (I) = - jH [QfJ(^z;a) + ^«,5;cr)] F <'> (a|(z 2 - 0"^) d£, (3.52) 

where x = r(^ ) 2 /(9Ag)i , 

<?<">(*„*») = - C, + 1 l«2) -*(?)]}, (3.53) 

F (/) is given by (A 17), C\ is given by (A 7) and xp denotes the digamma function defined 
on page 258 of Abramowitz & Stegun (1964). 

The pressure distribution P D and the particle displacement A can now be found from 
the coupled equations (3.29) and (3.51). These equations, together with (3.15), (3.16) and 
(3.20), imply that 

U D ~ax 2 f"Hz and V D ~-df'Hz as z 2 — + 0, (3.54) 

while 

P D _ „ j ln ° ~ C i + 3 fyl 2 ) - ^( 3 )] | as -» 0 (3.55) 
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which shows that the expansions (3.7)-(3.10) become invalid (as already anticipated) when 
x 2 becomes sufficiently small, i.e. when 

X = (TX = 0(1). (3.56) 

The relevant solution is constructed in appendix B where it is shown that the pressure 
fluctuation that matches onto (3.55) as x — + oo remains bounded at x = 0. 

Equations (3.20), (3.29) and (3.51) also imply that 

A ~ - ~ ^ 2 — |ln a — Ci| as x 2 — *■ oo, (3.57) 

which, when combined with (3.15) and (3.45), shows that the viscous wall-layer thickness 
increases like x% as x 2 — > ► oo and that it fills the entire Blasius boundary layer when 
xJRh = a 3 x 2 = 0(1). However, (3.7) implies that the streamwise velocity component of 
the vortex flow will still be small compared to the Blasius profile in this region. This, in 
turn, shows that the vortices undergo their entire life cycle from initial algebraic growth to 
ultimate viscous decay while still remaining small compared to the basic Blasius flow. 

4. Structure of the instability waves 

We now consider the oblique instability modes that are generated by the excitation 
device and which initially grow in accordance with linear dynamics. If it were not for the 
streamwise vortices, only the slowly growing Tollmien-Schlichting waves would be amplified 
by the mean flow but, as noted by Prandtl (1935) and Stuart (1965), the rapidly growing 
streamwise velocity perturbation can cause the mean flow to become inflectional and thereby 
support inviscid Rayleigh instabilities that exhibit much larger growth rates (when a is 
sufficiently large relative to .ft -1 ). The inflection point will always lie close to the wall 
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since the vortex velocity is always small compared to t/jg (the actual distance turns out 
to be 0(a)). The maximum inviscid growth of the instability waves will take place in the 
streamwise region where x-i = 0(1) since it follows from (3.44), (3.50), (3.52) and (3.57) 
that 

Uyy — 1 ► 0 as X2 —>■ oc at Y = constant. (4-1) 

The expansions (3.7)-(3.9) and (3.33)-(3.35) and the inequality (3.2) show that the 
cross-flow velocity components of the vortex remain small compared to its streamwise com- 
ponent in this region. This means that the initially linear instability modes, whose wave- 
lengths are small on the x 2 scale, are only affected by this latter velocity component and 
therefore behave (to lowest order of approximation) like linear perturbations of the uni- 
directional transversely sheared mean flow U(xi, y, z; a), with the slow streamwise variable 
X 2 entering only parametrically, i.e. playing the role of a constant. It therefore follows from 
inviscid instability theory that the wavelength of the most rapidly growing waves will be 
0(Sa~ 1 ) which is long compared to the boundary-layer thickness but small compared to 
the downstream distance from the excitation device, say x = a 3 R^C where l is order one. 

The initially linear non-interacting instability waves will eventually become nonlinear 
when their amplitudes become sufficiently large. However, the resulting nonlinear effects 
will be confined to the critical layer since the linear growth rates are small compared to the 
corresponding wavenumbers and the Reynolds number is large. 

The flow outside this layer will still be determined by linear dynamics. Its velocity and 
pressure will expand like 

u = U + e(u, av, !»)+•••, (4-2) 
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p = P + cp + • • • , 


(4.3) 


where u, v, w, and p are functions of 

x = ax - a 4 R$t, (4.4) 

y, z and 

t = at , (4.5) 

and e< 1 characterizes the local amplitude of the unsteady disturbance in the streamwise 
region where nonlinear effects first become important. The precise relationship between e 
and a will be specified below when the flow in the critical layer is considered. 

Substituting (4.2) and (4.3) into (2.3) and (2.4) and linearizing the result about the 
mean flow U yields 

D(u, v,w) + (U y v + U s w,Q,0) + {Px,°~ 2 Py,Pz) = 0{o~ l R~ i), (4.6) 

+ Vy + Wz = 0, (4.7) 

where D = d/dt + Ud/dx is the leading-order convective derivative relative to the mean 
flow. These equations are just the familiar equations for the linear perturbations about a 
uni- directional transversely sheared mean flow (Goldstein 1976; Henningson 1987). It is 
well known that the velocity fluctuations can be eliminated between them (see Goldstein 
1976, pp. 6-10 for a detailed derivation) to obtain the following equation for the pressure 
fluctuation 

D (<r~ 2 pj, y + Vlp) - <T~ 2 2U y Pxy - WzPxz = 0(a" 1 f?"i), (4.8) 

where = d 2 jdx 2 + d 2 jdz 2 is the Laplacian in the horizontal plane. 
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Since a < 1, our interest is in the long-wavelength instability wave solutions to (4.8). 
We subsequently show that, as in the case of a strictly two-dimensional mean flow, the rele- 
vant x-scale growth rates are 0(a 4 ), which means that the solutions grow on the relatively 
long scale 

Xj = a 4 (x - a 3 R^t) = c 3 x. (4.9) 

As a minimum, we want a to be large enough so that these growth rates are at least as 
large as the relatively small, i.e. 0(R~ ^), viscous growth rates corresponding to the upper- 
branch scaling for the Tollmien-Schlichting waves (Bodonyi & Smith 1981; Goldstein & 
Durbin 1986). This means that we should require 

R~’k=0{cr) (4.10) 

which is consistent with the inequality (3.2) and includes the case R~& = o(a) in which the 
growth rates are larger than the corresponding upper-branch growth rates of the Tollmien- 
Schlichting waves. In any event, the 0(cr~ 1 R ~ $) error terms will then be small compared 
to the wave-growth terms in (4.6)-(4.8) and these equations will then be accurate enough 
to determine the instability wave solutions outside the critical layer - provided we allow for 
a thin Stokes layer in the near wall region where y = 0(a~ 1 R~%). 

As already indicated, the initial upstream disturbance is best represented by a pair 
of equal- amplitude oblique instability modes having the same streamwise wavenumber and 
scaled frequency but opposite spanwise wavenumbers. These modes form a standing wave 
in the spanwise direction that exhibits its most rapid growth when its spanwise wavelength 
is twice that of the mean-flow distortion, i.e. when its spanwise wavelength is 2n /(}. 

The resulting solution for the unsteady portion of the flow will then be of the form 
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originally assumed by Henningson (1987), namely 

( u,v,w,p ) = Re v,w,p)e' x + (4.11) 

where Al(a;x) is a slowly varying amplitude function that accounts for the growth of the 
instability waves, Re z)] is a spanwise- variable mean-flow distortion that is gen- 

erated by nonlinear effects in the critical layer, 

X = a(x - act), (4.12) 

and the real quantities a(cr) and c(a) are the scaled streamwise wavenumber and phase 
speed correct up to but not including 0(a 3 ) terms, a and c possess expansions of the form 

a = Uo b • • ' ) c = Cq + • • • (4.13) 


as a —* 0 where ao and c 0 are order-one constants. The corresponding scaled Strouhal 
number or (non-dimensional) angular frequency is 

S = ac (4.14) 


where S is an order-one real constant. Substituting (4.11) into (4.8) shows that, outside the 
Stokes layer, the function p of zi, y and z is determined to the required order of accuracy 

by 


Py 1 

I 

0 2 pz 

(U-c) 2 , 

I 

y 

[(U-c) 2 


a 2 a 2 p 

(U-c) 2 


= 0 


(4.15) 


where 

q = a - a 3 iA^'/A^ and c = aS/a, (4-16) 

which insures that the instability wave exhibits only spatial growth. It follows from (4.6), 
(4.7) and (4.11) that the velocity fluctuations are determined in terms of p by 

UyPy U- Z p- Z p 


u = — 


a 2 a 2 (U - c) 2 a 2 (U - c) 2 U-c ’ 


(4.17) 
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(4.18) 


v = 


w — 


1 Py 

o 2 a(U — c)’ 
ip* 


a(U — c) 

The solution to (4.15) that satisfies the homogeneous boundary condition 


(4.19) 


p — > 0 as y — ► oo 


(4.20) 


and matches onto the Stokes-layer solution as y — ► 0 is derived in the following section. 


5. Solution outside the critical layer 

5.1. The Stokes layer 

Viscous effects, if they enter the reckoning, can only affect the solution within the critical 
layer and within a thin Stokes layer induced by the no-slip boundary condition at the wall. 
The appropriate scaled transverse variable for the latter region is 

Y = oR*y. (5.1) 

Substituting (4.2), (4.11) and (5.1) into (2.3) and (2.4) shows that the unsteady flow is 
determined to the required order of approximation in this region by 

— i oS(u, oR^v, w ) + (iap, R^py,pz) — <7(u, oR^ v, w)yy, (5.2) 

idu + crRivy + = 0, (5.3) 

together with the no-slip boundary condition 

u = v = w = 0 at f = 0. (5.4) 

The velocity fluctuations can be eliminated between (5.2) and (5.3) to arrive at an equation 
for the pressure fluctuation. The equation for the transverse velocity fluctuation v can then 


20 



be obtained by combining the Y derivative of this equation with the transverse component 
of (5.2). The solution that satisfies (5.4) and does not exhibit exponential growth as Y — * oo 
is 

v = a~ x RT^ au~ 3 (u>Y + e~ w ^ - 1 ^ (5.5) 

where u> = e~‘T and a is an arbitrary function of xj and z. It follows from (5.5) and the 
transverse component of (5.2) that 

Py = aR-$a(Y -u> -1 ) (5.6) 

which now can be used to derive the appropriate boundary condition for (4.15) as y — ► 0. 

5.2. The inviscid wall layer 

Introducing the wall-layer variable (3.32) into (4.15) and (5.6) and integrating the former 
with respect to Y subject to matching with the latter as Y — ► 0 shows that 

p = ab -f 0(a 5 ) (5.7) 

for y = 0 ( 1 ) where b is an arbitrary function of x\ and z that has an expansion in a 
containing terms up to but not including 0(a 4 ). It turns out that 

b = bo + • • • + icr 3 i> 3 ,- + — (5.8) 

as a — ► 0 where 6 o depends on z only and the coefficients 6 0 to 63 , are purely real so that 
60 and <J 3 bzi are the leading-order approximations to the real and imaginary parts of 6 . 

To determine the dispersion relation, the small-cr expansion of py is needed. By using 
(5.7) together with the expansions (3.33) and (4.16), (4.15) can be integrated with respect 
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to Y to obtain, after some manipulation, 


PV 

( U — CTc) 2 


= cr 3 P 


+ i V 

A o 


A 2 y c (y - y e ) J 

^ + p.( < l> ± -<t>~ ) 


cY 2 iA*' 


+ 


*8 


+ 


«y c 2 (y - y c ) 2 At T \ 

where V = b s d/dz + b ig — a 2 b, the transverse position of the critical level y c 
determined by the condition 

U{i,y c ,z-,c) = ac, 


(5.9) 
<ry c is 

(5.10) 


the ± superscript indicates different values for Y ^ Y c , <f> is the real function of Y and z 
determined by 


\{Y - Y c ) 3 $y = 0(£,Y,z;c 7 ) - U(t,Y c ,z;a) - ±X &Y 4 - Y‘ «) - \X 0 C(Y - Y c ), (5.11) 
with ^ = 0 at y = 0, and 

P= U YY {t,Y c ,z-a)-\\lY? (5.12) 

is the normal derivative of the scaled mean vorticity at the critical level. Matching (5.9) 
with the Stokes-layer solution (5.6) as Y — > 0 shows that 


a = -bn + a 2 b 


(5.13) 


where the small-cr expansion 

Y c = c/X 0 + 0(cr 3 ), (5.14) 

which is easily obtained from (3.7), (3.11) and (5.10), has been used. 

It follows from (4.17) and (5.7) that the discontinuity in (5.9) results in a jump in the 
streamwise component of the velocity fluctuation 


A u - —a 3 V 


ft(<P + ~ 4> ) /X la 1 1 


+ 


(5.15) 
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across the critical layer. Matching this jump with the one induced by the flow in the critical 
layer determines the integration ‘constants’ <f>^ (which are at most functions of x\ and z ). 
The velocity jump corresponds to a logarithmic phase shift of n when the critical layer is 
linear, which leads to the requirement that 

p + — cf)~ — + i 7r as x\ — ► — oo. (5.16) 

This ensures that the nonlinear solution (to be discussed below) will match onto the appro- 
priate linear solution in the upstream region. 

5.3. The main boundary layer 

Equations (5.7) and (5.8) together with the expression for p y obtained by substituting 
Y = y/a into (5.9) and re-expanding the result suggest that, in the main part of the 
boundary layer, p should expand like 

p = a (b 0 + ap 2 + • • •) + i<r 4 (63, ■ + <rp 5 i + • • •) (5.17) 

where p 2 and ps, are real functions of a:i, y and z and only the first two terms in the small-a 
expansions of the real and imaginary parts of p are shown. Substituting (3.11) into (3.7), 
inserting the result together with (5.17) into (4.15), and integrating with respect to y yields 

P2y = d 2 f'\ p5iy = d hi f'\ (5.18) 

where d 2 and d$i are at most functions of x% and z. Matching with the real and imaginary 
parts of the wall-layer solution (5.9) then shows that 

d 2 = V o(-r4-), (5.19) 
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(5.20) 


^ G^) + T ° 


H . ■ . \A^' 1 

Al ( ^ ~ 4> } " A^4t + a * R \<fo 


0"J 


where 


Vo — bo-d/ dz + bo ii — alb 0 , 


(5.21) 


and 


V 3 i = b 3ii djdz + b 3iii - alb 3i + 2a 0 Im(i,4 t 7>l t )&a, 


(5.22) 


are the scaled leading-order approximations to the real and imaginary parts of V. 


5.4. The outer layer 

Equation (5.18) shows that the solution (5.17) does not satisfy the free-stream boundary 
condition (4.20) and it is, therefore, necessary to construct the solution for p in the outer 
region described by (3.21). Substituting (3.21) into (5.17) and making use of (3.12) and 
(5.18) shows that, in this region, p expands like 

p = a pi + ■ • • + i <J 4 p4i H (5.23) 

where p\ and p 4 , are real functions of xj, y and z and only the leading-order terms in the 

small-<7 expansions of the real and imaginary parts of p are shown, pi and p 4t - must satisfy 

Pi = b 0 , pi s = d 2 , p 4 i = b 3i , p Ai - = d 5i at y = 0 (5.24) 

in order to match with the main-boundary-layer solution (5.17) and (5.18) and 

Pi, pAi ~ ► 0 as y oo (5.25) 

in order to satisfy the free-stream boundary condition (4.20). Equation (3.22) shows that 
the mean streamwise velocity U approaches unity at a fast enough rate to ensure that (4.15) 
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reduces to the Helmholtz equation when y = 0(1). Therefore, substituting (4.16) and (5.23) 
into (4.15) leads to 

Plyy + Plzz ~ «0Pl = 0, (5.26) 

Piiyy + pAizz - «oP4.- = -2a 0 Im(i J 4 t/ /A t )pi. (5.27) 

Up to this point, no restrictions have been placed on the z dependence of the solution. 
However, as indicated above, our interest is in the case where the upstream linear solution 
has a standing wave behavior in the spanwise direction with a spanwise wavelength twice 
that of the imposed cross-flow velocity Wo. This means that b should be of the form 

b — 2a cos 0z + 0(cr 3 ) (5.28) 

where the normalization 2a has been inserted for convenience. 

It now follows from (5.19) and (5.24)-(5.26) that the relevant solution to (5.26) is 

pi = 2aoe - ^ cos fiz (5.29) 

where 

fc= (al + P 2 )* = W (5.30) 

The boundary-value problem (5.24), (5.25) and (5.27) only possesses solutions for certain 
values of Im(i A^' /A^) since pi is a homogeneous solution to (5.27). These values can be 
found without explicitly solving for by integrating the difference between pi times (5.27) 
and p 4 i times (5.26) from y = 0 to oo, integrating the result from z = 0 to 2n/fl and then 
using (5.24), (5.25) and the z-periodicity of p and U to arrive at the following solvability 
condition 

r2njp ftn/P roo 

/ (6(^5. - d 2 b 3 i)dz = 2a 0 Im(iA t /A t ) / / pjdydz. (5.31) 

Jo Jo Jo 
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It now follows from (5.19)-(5.22), (5.29) and (5.30) that 

( cos 6 4 - (A^ - i«,A t ) = 

\ cose?/ 

Q r2n/fi =2 _ \2 

-i- 7 — / -^-p. (<f> + - <j>~) (l + cos20 cos2fiz) dz + -. .P A^ (5.32) 

2nJ 0 Xq v ’ <r 5 Rl(2S)l 

where 9 = arctan(/3/d 0 ) is the obliqueness angle of the instability wave and the real constant 
k, has been introduced to account for the 0 (<r 3 ) term in the expansion of the streamwise 
wavenumber in the upstream linear region. is completely determined by the linear prob- 
lem outside the critical layer but its explicit form is not given here because it is not needed 
for the subsequent analysis. 

Equations (3.44), (3.45), (5.12) (5.14), (5.16) and (5.32) imply that 

A * — ► Age** 1 as zi — ► -00 (5.33) 


where Aq is a constant, k = n r -f i/c,-, 


COS 9 n C„ - „ „ 7T ci A? 

K r = -—2 Q cos 20 77 ® + . - ° _- r 

cos 2 0 + l A 0 4A 0 a 5 Rl(2S)% 


is the common initial parametric growth rate of the oblique modes and 


3 r 2 ”/^ - 
n = — Uyy(£, Y c , 2 ; a) cos(2/3z)di 

In Jq 


(5.34) 


(5.35) 


with U given by (3.44), (3.50) and (3.52). The last term in (5.34) accounts for the viscous- 
Stokes-layer effect and is negligible when a > R~& . The first term in (5.34) represents 
the growth produced by the resonant interaction with the streamwise vortices. It is always 
possible to make k t positive by appropriately selecting the imposed cross-flow velocity Wq 
which is now assumed to be the case. 
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6. Nonlinear-critical-layer effects 

Since (5.33) and (5.34) imply that the linear instability wave continues to grow as it 
propagates downstream, nonlinear effects eventually come into play and, as already noted, 
this first occurs locally within the so-called critical layer. The thickness of the critical layer 
(on the y scale) turns out to be of the same order as the growth rate so the appropriate 
scaled transverse coordinate for this region is 

fj = (y - crc/A 0 )/<r 4 = (Y - c/X 0 )/a 3 . (6.1) 

Nonlinear terms produce a critical-layer velocity jump at the same order as the linear/ 
parametric-growth terms when the amplitude scale c, which was defined in (4.2), is chosen 
to be 

€ = < 7 10 . ( 6 . 2 ) 

Viscous effects will enter into the dominant balance for the critical layer while making only 
insignificant modifications to the external flow when the Benney-Bergeron parameter 

A = l/a 13 fii (6.3) 

(Benney & Bergeron 1969) is order one. The implied wavelength-Reynolds-number scaling 
(6.3) is more restrictive than (4.10) in the sense that the Stokes-layer contribution to (5.34) 
is negligible when A = 0(1). However, retaining this term in (5.34) while assuming A = 0(1) 
in the critical-layer analysis leads to a kind of distinguished scaling from which the more 
viscous case (4.10) can be recovered as a limit (see below) and we therefore adopt (6.3) as 
the appropriate wavelength-Reynolds-number scaling for the present problem. 
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Equations (5.10) and (6.1) imply that 


Y - Y c = a\f, -%) + ■" (6.4) 

where 

Vc(z\ a) = -Aq 1 U(i, Y c , z ; a) + ^A 0 F C 4 + \IY C (6.5) 

accounts for spanwise variations in the position of the critical level. Introducing (6.4) into 
the expressions for it, i), w and p obtained from (4.17)-(4.19), (5.7) and (5.9), re-expanding 
the result, and using (3.33)-(3.35), (4.2), (4.3), (4.11), (6.2) and (6.3) shows that the critical- 


layer flow should expand like 

u = ac + (t 4 \ 0 (v - fj c ) + cr 7 ui + ct 10 u 2 + <t 13 u 3 H , (6.6) 

v = a 12 v i + a 15 v 2 + a 18 v 3 + • • • , (6.7) 

w = o 7 w\ + a 10 w 2 -f <t 13 uj 3 + • • • , ( 6 . 8 ) 

p — a~ 2 R~ x Po + <7 n 2a cos/3zRe(AV*) -f a 14 p 2 + <z l ‘ps + • • •, (6-9) 


where the u n , v n , w n , and p n are functions of X , x \ , t) and z that have expansions in a 
containing terms up to but not including 0(cr 3 ). Matching (6.6)-(6.9) with the external 
linear solution requires that 

ui - [U Y (t, Y c , z- a) - ±\ 2 Y? - \\ 0 t\ f) 

+2f3 2 cos j3zRe(A^e ,x )/ \oa(fj - fj c ), (6.10) 

and 

wi — ► -2@ sin f3zRe(iA^e >x )/\o(f] - ijc), (6.11) 
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as 77 — * ±00. Matching the critical layer induced velocity jump with the external jump 
(5.15) requires 


1 r + 00 r 2 n . 1 

-J J q e~ ,x u 3n dXd7j = -V 0 [/i(0 + - <A~)/A^] A f . (6.12) 


Multiplying (6.12) by cos j3z, integrating from z = 0 to 2 tt //3 and then combining the result 

with (5.32) yields the jump condition 

8 r 2n /P /+°° r 27r -y 

77— j / e ,x u 3 -cos8zdXdr]dz = 

Jo J — 00 Jo 

h 


. Ao 
1— 
«o 


( C0S 9 + cra») ( '* t ' " iK ' " A “ (l?) 


(6.13) 


Notice that, for the reasons given above, the asymptotically small Stokes-layer term has 
been retained in this result. 

Substituting the expansions (6.6)-(6.9) into (2.3) and (2.4) and combining the result 
with (6.10), (6.11) and (6.13) leads to a boundary-value problem for the u n , v n , w n , and 
p n . When this is expressed in terms of the strained coordinate 


V = V ~ ic 


(6.14) 


and the corresponding velocity component 

v n = v n - fj c -w n , for n = 1,2,3, (6.15) 

we find that 

Vi = -2(d 2 + (3 2 )cos(3zRe(iA^e i;f )/Ao, (6.16) 

and, more generally, that the entire problem is just the viscous version of the one solved in 
Goldstein & Choi (1989), a special case of the one considered in Goldstein & Lee (1992), and, 
except for differences in notation, precisely the one considered in Wu, Lee & Cowley (1993). 
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The spanwise variation of the mean flow (which is the new feature in the present analysis) 
now enters the problem only through the transverse boundary condition (6.10) as well as the 
transverse boundary conditions for ui, U 3 , V 2 ^3 and W 3 . But, since the spanwise variation 
of the mean flow introduced by these conditions only affects the linear/parametric-growth 
terms, the solution to the present critical-layer problem can be easily deduced from the 
solutions given in the previous investigations. Substituting this solution into the velocity- 
jump condition (6.13) leads to the following equation for A' 

A*' = KA' + \M r /*’ /$'(®il6,6M t (6M t (6M r (6 + 6-xi)d6d6 (6.17) 

J — CO J — CO 

where 

M = 27rAQ<5o cos 2 0/(cos 2 8 + 1), (6.18) 

the asterisk denotes complex conjugation and the initial parametric growth rate n r is given 
by (5.34). The kernel function K is fairly complicated in the general viscous (i.e. order-one 
A) case considered by Wu, Lee & Cowley (1993) and Leib & Lee (1994), but, in the inviscid 
limit first considered by Goldstein & Choi (1989), it is a relatively simple polynomial of the 
streamwise coordinate and is given by 

K - \ tan 2 0cos20(zi - fr) [cos20(:ri - fi )(6 - 6) - (*i - fi ) 2 - (*i - 6) 2 ] • (6.19) 

As already indicated the Stokes-layer contribution to k t is negligible when A = 0(1). How- 
ever, retaining this contribution and taking the limit as A/k^ -+ oo of the integral term in 
(6.17) as was done by Wu, Lee & Cowley (1993), leads to the amplitude-evolution equation 
corresponding to the upper- branch Blasius-boundary-layer scaling a R& = 0(1). In the 
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present notation, this equation is 


A^' = kA^ — iA $Mta,n 2 9s\n 2 6cos29(jg)$T(^)A^ f |At(£)| 2 d£, (6.20) 

J — OO 

where A = XqOIqX/c^ is a rescaled Benney-Bergeron parameter. 

7. Results and discussion 

In one sense, the oblique-mode amplitude equation (6.17) (or its highly viscous limit 
(6.20)) can be thought of as the final result of this paper. While this equation is now well 
known, its application to the present situation is new as is the formula (5.34) for the initial 
parametric growth rate of the oblique modes k t . The latter accounts for the parametric- 
resonance effects that allow the oblique modes to grow faster than the two-dimensional 
waves which, in turn, allows the oblique-mode interaction described by (6.17) to become 
the dominant interaction in the initial nonlinear stage of the transition process. The initial 
parametric-resonant interaction takes place between the relatively weak spanwise-periodic 
mean-flow distortion (i.e. the streamwise vortices) and a pair of equal-amplitude oblique 
modes that form a standing wave in the spanwise direction with spanwise wavelength equal 
to twice that of the mean-flow distortion. Equation (5.34) shows that its spatial growth 
rate is enhanced by the inflectional nature of the mean velocity profile. 

The transverse position of the mean-flow inflection point y, = oY a is determined by 

Uyy - \*lY 2 =0 at Y = Y s . (7.1) 

Figure 1 shows how Y, varies with x<i for an initial spanwise velocity that behaves like 

B(z) = W o (0, z ) = (2$)~ 1 Bq sin 2 0z (7.2) 


31 





at the wall. For this class of imposed cross flows, (3.20), (3.29) and (3.51) can be solved 
analytically to obtain 

A = Z (-X)*- 1 r Q { r } (Z, -z\ <T)e-^-«df , (7.3) 

m=4 J0 

and 

Pn- Z - Z = W X2X2 , (7.4) 

where x = 9AQ/[2/?r(i)] 3 . The 0(l/ln£x) terms in (3.44) and (3.53) were kept when 
computing the results presented here since these terms can be significant even at relatively 
small values of a. The mean-flow distortion U initially grows linearly with increasing x 2 (see 
(3.7), (3.15), (3.33), (3.44), (3.52) and (3.55)) causing Y s to move out from the wall. This 
type of linear disturbance growth is similar to the algebraic growth proposed by Ellingsen 
& Palm (1975), Hultgren & Gustavsson (1981) and Landahl (1990) as an alternative or 
‘bypass’ transition mechanism. However, only the mean-flow distortion (or vortex flow) 
undergoes this type of growth in the present study and the follow-on ‘secondary’ instability 
modes exhibit the more conventional exponential-type growth. In the present analysis, the 
initial mean-flow growth is eventually reversed by viscous effects once the viscous wall layer 
expands to fill the entire Blasius boundary layer, i.e. once x 2 becomes O(o~ 3 ). However, 
figures 1 and 2 show that the mean-flow inflection point produced by the distortion actually 
vanishes before this stage is reached. 

The corresponding initial parametric growth rate k t is shown in figure 3. Near the exci- 
tation device, i.e. at small values of (, the mean-flow inflection point lies very close to the wall 
and, as shown by the figure, K r is negative when <rR& = 0(1) and given by the Tollmien- 
Schlichting-wave result when aR& = 0(1). The actual unsteady flow in this region would 
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Figure 2: Scaled streamwise velocity profiles of the vortex flow within the viscous wall layer 
for a — 0.05, (3 = n /2A and various values of x^. 


34 
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Figure 3: Initial parametric growth rate vs. scaled Strouhal number for a — 0.05 and 
various values of l. (a) Bq — — 12 7t, (3 = tt / 24, A = 0; (b) B 0 = - ft = 7t/24, A = 0; (c) 
B 0 = —16 7r, /3 = ti/ 18, A = 0; (d) B 0 = f3 = ti/ 24, A = 1/a 3 . 
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probably be dominated by two-dimensional viscous instabilities, i.e. by Tollmien-Schlichting 
waves, with the three-dimensional inviscid instabilities emerging further downstream where 
the mean-flow distortion has become sufficiently large. The continued growth of this dis- 
tortion can then lead to much larger growth rates for the three-dimensional instabilities 
than those of the Tollmien-Schlichting waves. This is consistent with the experimental ob- 
servations, which probably correspond to the case where crR^ = 0(1), and show that the 
three-dimensional structures (which can be identified with the oblique modes described by 
the present analysis) gradually emerge from the initial two-dimensional motion. Figures 3 
and 4 suggest that these structures will eventually decay once the mean-flow inflection point 
has been eliminated by viscous effects. However, the oblique modes will most likely become 
nonlinear before this occurs and their amplitudes will then determined by (6.17). 

Since the nonlinear critical-layer equations (when expressed in terms of f j) explicitly 
involve the spanwise variation of the critical-level position, it would have been impossible to 
anticipate that the nonlinear oblique-mode amplitude would be determined by (6.17) which 
had previously been derived only for nonlinear interactions on strictly two-dimensional mean 
flows. However, the transformation (6.14) and (6.15) eliminates this spanwise variation from 
the problem and the two-dimensional mean-flow equation (6.17) is therefore obtained. This 
equation can be rescaled to eliminate the parameters k and M (see Goldstein & Choi 1989; 
Wu, Lee & Cowley 1993; Leib & Lee 1994) so that the resulting solutions depend only on 
the obliqueness angle 9 and the scaled viscous/growth-rate parameter X/nf. 

The numerical results for the rescaled oblique-mode amplitude are relatively universal 
and effectively the same as those given by Wu, Lee & Cowley (1993) in the viscous case 
and by Goldstein & Choi (1989) in the inviscid limit. We therefore do not present any new 
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computations for A\ but merely note that, as was shown by Goldstein & Choi (1989), the 
solutionsN to (6.17) become singular at some finite downstream position, say x„, and that 
the local asymptotic expansion in the vicinity of this singularity is 

~ a,(x. - zi)- 3 ~ iv (7.5) 

where a„ and <p are constants that aTe determined explicitly in Goldstein & Choi (1989). 
This suggests that the initially linear instability waves can exhibit nonlinear break down 
before the mean-flow distortion (which supports these waves) is able to decay. This is 
analogous to the argument used by Cowley (1987) and Wu, Lee & Cowley (1993) to explain 
the breakdown of Stokes layers which (unlike the experimentally observed flows) would 
exhibit global stability over each oscillation cycle if nonlinear breakdown did not occur. 

It is easy to solve the highly viscous amplitude equation (6.20) analytically and thereby 
show that the oblique modes continue to grow linearly in this limit. However, this solution 
shows that the wavenumber correction Re(<7 4 j4^/i>f^) increases exponentially with increas- 
ing Xi - indicating that the assumed wavenumber scaling must eventually break down. The 
next stage of evolution should be governed by the full non-equilibrium equation (6.17). The 
highly viscous solutions to this equation indicate that the explosive growth still occurs in 
this case (Dr. Sang Soo Lee, private communication) - which suggests that explosive growth 
will occur even when aR& = 0(1). 

The amplitude B 1 of the spanwise- variable mean-flow distortion Re[B^(xi)g(y, £)] in- 
duced on the external linear flow by the nonlinear critical-layer interactions is given by 
Goldstein (1994) as 

fit = b 3 r (*! - 6)e-^>- 6)3 U f (6)| 2 dfcdfc, (7.6) 

J — OO J —oo * 1 
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which shows that B t also becomes singular when x\ — t - x s and behaves like 

^ t ~ro*»|a.| 2 (*.-*ir 3 (7-7) 

in this limit. 

Equations (3.7), (3.15) (3.19), (3.33), (4.2), (4.11) and (6.2) show that this induced 
spanwise-variable mean-flow velocity exceeds the initially imposed distortion velocity (in 
both the wall layer and main boundary layer) when 

x a - x\ < ct 2 . (7.8) 

The corresponding large amplitudes can occur without violating the present asymptotic 
scaling because the instability wave growth is self-induced and does not depend on the 
initial parametric growth once the oblique-mode amplitude becomes sufficiently large. 

Of course, the oblique-mode amplitude and associated spanwise-variable mean-flow dis- 
tortion cannot continue to increase indefinitely, and a new stage of evolution must eventually 
be reached when, as pointed out by Goldstein & Lee (1992), the growth rate Re(a 4 A^/^) 
becomes of the order of the wavelength scale <7, i.e. when 

x s - xi = 0(<r 3 ), (7.9) 

which is much shorter than the length scale (7.8) at which the induced mean-flow distortion 
exceeds the one imposed on the flow. The unsteady flow and spanwise distortions are now 
0{oUoo ) in this smaller region and evolve on the relatively short wavelength scale x. The 
resulting motion is, as noted in Goldstein & Lee (1992), primarily inviscid and determined 
by the triple-deck equations (Stewartson 1969; Messiter 1970), but with no viscous terms 
appearing in the nonlinear lower-deck equations. This latter region has the same thickness 
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as the wall layer associated with the original vortex system. There must, of course, be a thin 
viscous region underlying the triple-deck structure. The flow in this lower region, which 
is governed by the usual three-dimensional boundary-layer equations with the externally 
imposed pressure gradients determined by the flow in the triple deck, can, as pointed out 
by Smith & Burggraf (1985), undergo large scale separation. However, it would not be 
possible to determine the flow in the triple deck, even if separation did not occur, because 
the downstream boundary conditions for this effectively elliptic problem depend on the 
downstream flow which is, in essence, turbulent in the situation of interest heTe. 

While the relationship between the spanwise-wavelength scale a and the Reynolds 
number R has not been completely fixed in the present analysis, we have required that 

/Rl. (7.10) 

The mean-flow distortion would be governed by the full three-dimensional triple-deck equa- 
tions in the limit a — ► l/Rt. The present scaling corresponds to the short- wavelength limit 
of the triple-deck problem which was used by Rozhko & Ruban (1987) and subsequently by 
Choudhari, Hall & Streett (1992) to study Gortler vortices. However, the most appropri- 
ate distinguished scaling for the present problem corresponds to the order- one- wavelength 
limit a —* 1. We chose to consider the long-wavelength limit cr < 1 in the hope that the 
analytical solutions that were obtainable in this limit would lead to a better understanding 
of the physical mechanisms involved. 

Large scale numerical computations are required when a = 1 since the spanwise- 
variable mean flow (3.3) and (3.4) can only be determined by solving the parabolized 
Navier-Stokes equations (3.5) and (3.6). The relevant solutions still grow linearly with 
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increasing downstream distance when xi is sufficiently small. In fact they are still given by 
(3.7)-(3.10), but with a and a set equal to unity and with U D , V D , W D and P D given by 
(3.54), (3.17) and (3.55), respectively. This shows that the algebraic growth of the steady 
distortion is still important when a — 1. 

While the initial linear instability wave can still be determined from the generalized 
Rayleigh equation (4.15) once the mean-flow velocity U is known, the relevant solutions must 
now be found numerically since a = 1. However, the long-wavelength solutions of (4.15) 
show that there is an effective upper branch to the neutral stability curve (see figure 3) - 
primarily because the mean-flow distortion ultimately decays to zero and thereby causes the 
inviscid instability wave growth rate (which is produced by this distortion) to vanish long 
before this decay is complete. Since this behavior also occurs when a = 1, all instability 
waves with sufficiently small initial amplitudes will now remain linear until their growth 
rates become sufficiently small. (Recall that the nonlinear amplitude of the instability waves 
scales with the growth rate and nonlinearity therefore occurs at smaller amplitudes when 
the growth rates become smaller.) 

The nonlinearity will then occur within a thin critical layer with the solutions outside 
this layer still given by (4.11) and (4.15)-(4.20). Hall & Horseman (1991), Horseman (1991) 
and Hall & Smith (1991) have studied the local critical- level behavior of the solutions to 
the generalized Rayleigh equation and have shown that it is nearly identical to that of the 
usual three-dimensional solutions to the Rayleigh equation for strictly two-dimensional mean 
flows. In a sense, all solutions of the Rayleigh problem for spanwise- variable mean flows 
behave like three-dimensional solutions to the Rayleigh problem for strictly two-dimensional 
mean flows in the vicinity of their critical levels. This is most easily shown by expressing 
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(4.15) in orthogonal curvilinear coordinates with one set of coordinate surfaces correspond- 
ing to surfaces of constant mean-flow velocity U - as was done, for example, by Goldstein 
(1976, pp. 6-10). 

The similarity of the solutions within the critical layer itself is even more dramatic. The 
critical-layer nonlinearity is still weak in the sense that it enters through an inhomogeneous 
term in a higher-order problem rather than through a coefficient in the lowest-order or 
dominant-balance equation and the instability wave amplitude At(xi) can therefore still be 
determined from a single amplitude-evolution equation, which is again given by (6.17) but 
with the inviscid kernel function K (which is still a simple polynomial of the streamwise 
coordinate) given by a slight generalization of the result (6.19). 

Appendix A. Solutions for C/ (f) and W^ l) 

In this appendix, solutions for (7 (/) and W (,) are constructed by first substituting 
(3.44) into (3-37) and (3.39)-(3.43) and then taking the Laplace transform of the resulting 
equations with respect to x 2 to obtain 

\ 0 sYW^ + P Dz - = W&, (Al) 

A osYUP - A 0 W { s n = Uyyy , (A 2) 

W (I) = U (I) = Uyy = 0 at Y = 0, (A3) 

and 

(jm _> A o (0j 6) + -4) as F^oo, (A 4) 

where 

(')= fe- ,IJ (')< 1*2 (A 5) 

Jo 
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denotes the Laplace transform with respect to X 2 , 


Q {m \s,z-a) = aX 0 l s ^ |ln [ct(A 0 s) i] - C\ + 5^(2)} B(z), (A 6) 

C\ = lim (F — In tj), (A 7) 

rj— *oo v ' 

and xj> denotes the digamma function defined on page 258 of Abramowitz & Stegun (1964). 
The relevant solutions to (A 1) and (A 2) are 

W (,) = -n(X oS )~^P Di [Gi(r?) - 3 - ^Ai(f/)j , (A 8) 

and 

U (1) = X 0 n(X 0 s)-$P D -- Gi(f))-Z-tAi(f,)-tfT(\)^£ A\(t)dt , (A9) 

where f) = (Aos)^T and Ai and Gi are the Airy functions defined on pages 446 and 448 
of Abramowitz & Stegun (1964). It follows from (A4), (A 9) and the asymptotic behavior 
of the Airy functions for large values of their argument (Abramowitz & Stegun 1964, pp. 
448-450) that 

(9A^)-ir (\)P D . Z - Z = -$?> - sU (A 10) 

which can be inverted to give (3.51). 

Before inverting (A 8) and (A 9), it is convenient to rewrite these equations as 

w U) = -X^r(\)-\X 0 s)^ P d -&'\s,Y), (A 11) 

and 

U {1) = T{l)-\X 0 s)-^ P D --F^ n {s,Y), (A 12) 

where G (/) and F^ 1 are the Laplace transforms with respect to X2 of the functions of 
tj = (A0/X2 )^Y determined by 

G {1) " + \t] 2 G (,) ' = -1, + i^F' 7 *" + At/F^' = -\r,G U) ', (A 13) 
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with 


G U) _ fO) _ pi')" _ o at V = 0, (A 14) 

and 

G (I) -* 0, _♦ _3~Sr(i) 2 as r] — * oo. (A 15) 

It follows from these equations that 

G (/) = 3" »r(|) 7 ( r) + f r »e-*7(i -<)<*«. (A 16 ) 

Jo 

and 

f<'> = -a-Sixijrd)- 1 / T r *e-‘u(l, t)d* + G (/) (A 17) 

Jo 

where r = 7/ 3 /9 and U denotes the confluent hypergeometric function in the notation of 
Abramowitz & Stegun (1964, p. 504). Using (A 10) when inverting (A 12) then leads to 
(3.52). 

Appendix B. Mean-flow solution for x = O(cr) 

In this appendix, we briefly discuss the mean-flow solution in the streamwise region 
characterized by (3.56). The lowest-order main-deck velocity field is a simple re-expansion 
of the solution (3.15)-(3.17) but, as can be anticipated from (3.29) and (3.55), with the 
particle displacement A set equal to zero (i.e. it is now of higher order). Equation (3.55) 
also implies that the pressure in this region should expand like 

P = aAA~ l R~^ P D + ■ •• (Bl) 

where A = (T~i R~^ <Cl,A = l/lnA and P D is a function of x and z that has the implicit 
A dependence 

P D = P Q (x,z) + APi(x,z). (B2) 
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Equations (3.7), (3.15) and (3.56) imply that the streamwise velocity distortion is now 
0(R~ i), i.e. of the same order as the imposed cross-flow velocity that produced it. 

The main difference from the x = 0(cr 3 R$) solution comes in the viscous wall layer 
which now corresponds to 



Y = y/A = 0( 1). 

(B 3) 

The velocity in 

this region expands like 



U = A\qY + aA~ 1 R~^U , 

(B 4) 


V = (t&aa-'r-^v + ■■■, 

(B5) 


W = dA~ l R'^W H , 

(B6) 

where, like P D , 

the functions U, V, and W of x, Y , and z depend implicitly on A. 

The flow 

in this region is 

determined by 



\ 0 YUi+ \ 0 V + P D . = Uyy, 

(B 7) 


X 0 YWi + P Di = Wyy, 

(B 8) 


U £ + Vy+ W, = 0. 

(B9) 

with 

u = V = W = 0 at f = 0, 

(BIO) 


W -+ AB, U -* AxB’\n(AY) as Y oo, 

(B 11) 

and 

U — * 0, W -* AB as x — ► 0. 

(B 12) 
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The only difference from (3.36)-(3.41) is that the pressure gradient now appears in the 
streamwise momentum equation (B 7) and the particle displacement A does not appear in 
the boundary condition (B 11). This means that the pressure is completely determined by 
the solution in the lower deck and the now higher-order particle displacement is determined 
after the fact from the upper-deck problem. The solution procedure is essentially the same 
as before. Since the similarity variable tj remains order one when expressed in terms x and 
T, it is again appropriate to seek a solution of the form 

U = AxB'(z)F(ri) + U (I \x,Y,z-,A), W = AB(z)G(t]) + W^\x,Y,z\ A) (B 13) 

where the Laplace transforms of and are determined by (A 1) and (A 2) but with 
Y replaced by Y and with the boundary conditions now given by 

W (,) = U (f) = 0, = sP D - P D { 0, 5; A) at Y = 0, (B 14) 

and 

iy (/) ->0, U (I) -» \ 0 Qf\s,z-,A) as y -*■ oc, (B 15) 

where z\ A) determined from (A 6). 

It therefore follows that W (,) is still given by (A 8) but with fj = (Aos)^y, while U < - ,) 
is now 

U u) = A 0 7t(A 0 s)-^ { P D -- [Gi(77) - 3-i Ai(f))j 

-3*r(i) \P Dii + s 2 P d - sP D (0,z; A)] ± } . (B 16) 

Substituting this result into (B 15) and using the asymptotic behavior of the Airy functions 
for large values of their arguments (Abramowitz & Stegun 1964, pp. 448-450) shows that 

(9Ajj)-ir(i) [P Dgg + s 2 P D - sP d (0,z;A)} = -Q^\s,z-A), (B 17) 
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\ 


and since (B 7), (B 10) and (B 12) imply that P Di (0,z; A) = 0, this equation can be inverted 
to give 

(9Ag)-*r($ )(P Dii + P OH ) = A) (B 18 ) 

where Q^ m \x,z;A) is determined from (3.53). The solution to (B18) that matches onto 
(3.55) as x — » oo is 

rlZ j 0 r OO 

(9A5)-ir(I)P D = -£i 1 Q(%,C;4)Re{cot^ -< + i(*-01 

zn Jo Jo 

+ cotfi[z - C + i(i + ^)]}d^dC (B 19) 

where f denotes the Cauchy principal value. 
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